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Bimetric theory describes gravitational interactions in the presence of an extra spin-2 field. Pre¬ 
vious work has suggested that its cosmological solutions are generically plagued by instabilities. We 
show that by taking the Planck mass for the second metric, Mf, to be small, these instabilities can 
be pushed back to unobservably early times. In this limit, the theory approaches general relativity 
with an effective cosmological constant which is, remarkably, determined by the spin-2 interaction 
scale. This provides a late-time expansion history which is extremely close to ACDM, but with a 
technically-natural value for the cosmological constant. We find Mf should be no larger than the 
electroweak scale in order for cosmological perturbations to be stable by big-bang nucleosynthesis. 

We further show that in this limit the helicity-0 mode is no longer strongly-coupled at low energy 
scales. 


“The reports of my death have been greatly 
exaggerated.” 

—Metrics Twain 

Introduction. — The Standard Model of particle physics 
contains fields with spins 0, 1/2, and 1, describing mat¬ 
ter as well as the strong and electroweak forces. General 
relativity (GR) extends this to the gravitational inter¬ 
actions by introducing a massless spin-2 field. There is 
theoretical and observational motivation to seek physics 
beyond the Standard Model and GR. In particular, GR is 
nonrenormalizable and is associated with the cosmologi¬ 
cal constant, dark energy, and dark matter problems. To 
compound the puzzle, the GR-based A-cold dark matter 
(AGDM) model provides a very good fit to observational 
data, despite its theoretical problems. In order to be ob- 
servationally viable, any modified theory of gravity must 
be able to mimic GR over a wide range of distances. 

A natural possibility for extending the set of known 
classical field theories is to include additional spin-2 fields 
and interactions. While “massive” and “bimetric” theo¬ 
ries of gravity have a long history , nonlinear theories 
of interacting spin-2 fields were found, in general, to suf¬ 
fer from the Boulware-Deser (BD) ghost instability Q. 
Recently a particular bimetric theory (or bigravity) has 
been shown to avoid this ghost instability d, • This the¬ 
ory describes nonlinear interactions of the gravitational 
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metric with an additional spin-2 field. It is an extension 
of an earlier ghost-free theory of massive gravity (a mas- 
sive spin-2 field on a nondynamical flat background) [1-Q 
for which the absence of the B D ghos t at the nonlinear 
level was established in Refs. iffllll. 

Including spin-2 interactions modifies GR, inter alia, 
at large distances. Bimetric theory is therefore a can¬ 
didate to exp lain the accelerated expansion of the Uni¬ 
verse [H, [l^. Indeed, bigravity has been shown to pos¬ 
sess Friedmann-Lemaitre-Robertson-Walker (FLRW) so¬ 
lutions which can match observations of the cosmic ex- 
pansion history, even in the absence of vacuum energy 
[lj-l20jo Linear perturbations around these cosmologi- 
cal backgrounds have also been studied extensively |2^ 
[u| . The epoch of acceleration is set by the mass scale m 
of the spin-2 interactions. Unlike a small vacuum energy, 
m is protected from large quantum corrections due to an 
extra diffeomorphism symmetry that is recovered in the 
limit m —?► 0, just as fermion masses are protected by 
chiral symmetry in the Standard Model (see Ref. for 
an explicit analysis in the massive gravity setup). This 
makes interacting spin-2 helds especially attractive from 
a theoretical point of view. 

Cosmological solutions lie on one of two branches, 
called the finite and infinite branches 0 The infinite- 
branch models can have sensible backgrounds dill, 
but the perturbations have been found to contain ghosts 
in both the scalar and tensor sectors [H, [13, |4l| . Most vi¬ 
able background solutions lie on the hnite branch EMI. 


^ Stable FLRW solutions do not exist in massive gravity [2ll - [2^ . 

^ There is a third branch containing bouncing solutions, but these 
tend to have pathologies [ 4 III . 
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While these avoid the aforementioned ghosts, they con¬ 
tain a scalar instability at early times that in¬ 

validates the use of linear perturbation theory and could 
potentially rule these models out. For parameter values 
thought to be favored by data, this instability was found 
to be present until recent times (i.e., a similar time to the 
onset of cosmic acceleration) and thus seemed to spoil the 
predictivity of bimetric cosmology. 

In this Letter we study a physically well-motivated re¬ 
gion in the parameter space of bimetric theory that has 
been missed in earlier work due to a ubiquitous choice 
of parameter rescaling. We demonstrate how in this re¬ 
gion the instability problem in the finite branch can be 
resolved while the model still provides late-time acceler¬ 
ation in agreement with observations. 

Our search for viable bimetric cosmologies will be 
guided by the precise agreement of GR with data on all 
scales, which motivates us to study models of modified 
gravity which are close to their GR limit. Often this 
limit is dismayingly trivial; if a theory of modified grav¬ 
ity is meant to produce late-time self-acceleration in the 
absence of a cosmological constant degenerate with vac¬ 
uum energy, then we would expect that self-acceleration 
to disappear as the theory approaches GR. We will see, 
however, that there exists a GR limit of bigravity which 
retains its self-acceleration, leading to a GR-like universe 
with an effective cosmological constant produced purely 
by the spin-2 interactions. 

Bimetric gravity .— The ghost-free action for bigravity 
containing metrics g^i, and is given by @,0 


f A r Mi _ 

S = J d^x[-^^R{g)-^^R{f) 

+ Mp^AygVllC) + A/gJlm {g, ^i) ■ 


Here Mpi and Mf are the Planck masses for and , 
respectively, and we will frequently refer to their ratio, 


Mpi 


( 2 ) 


The potential R(X) is constructed from the elementary 
symmetric polynomials e„(X) of the eigenvalues of the 
matrix X = aJ g~^f, defined by 


X'^„X“,, = (3) 

and has the form [1,0| [3 

3 

V^y(X) = V5/3 o + Vs ^ /?„e„(X) + v7/34. (4) 


In the above, m is a mass scale and /?„ are dimension¬ 
less interaction parameters. /3o and parameterize the 
vacuum energies in the two sectors. Guided by the ab¬ 
sence of ghosts and the weak equivalence principle, we 
take the matter sector to be coupled to Then the 

vacuum-energy contributions from the matter sector 
are captured in Pq. We can interpret g^i, as the space- 
time metric used for measuring distance and time, while 
is an additional symmetric tensor that mixes nontriv- 
ially with gravity. As we discuss further below, the two 
metrics do not correspond to the spin-2 mass eigenstates 
but each contain both massive and massless components. 
Even before fitting to observational data, the parameters 
in the bimetric action are subject to several theoretical 
constraints. For instance, the squared mass of the mas¬ 
sive spin-2 field needs to be positive, it must not violate 
the Higuchi bound [HI, 0 , and ghost modes should be 
absent. 

In terms of the Einstein tensor, the equations of 
motion for the two metrics take the form 

Giiv{g) + w „ T^i,, (5) 

iWpi 

= 0 , ( 6 ) 

where are determined by varying the interaction 

potential, V. Taking the divergence of eq. (d) and using 
the Bianchi identity leads to the Bianchi constraint., 


= 0. (7) 

The analogous equation for carries no additional in¬ 
formation due to the general covariance of the action. 

Finally, note that the action o has a status similar 
to Proca theory on curved backgrounds. It is therefore 
expected to require an analogue of the Higgs mechanism, 
with new degrees of freedom, in order to have improved 
quantum behavior. The search for a ghost-free Higgs 
mechanism for gravity is still in progress [6ll| . 


The GR limit .— When bigravity is linearized around 
proportional backgrounds with constant c|3 


— 9n^ T 


Mpi 






ffAiA — c gp,u + 


( 8 ) 

( 9 ) 


the canonically-normalized perturbations can be diago¬ 
nalized into massless modes SG^iv and massive modes 

as [10 

5Gf_iy oc {5g^j,i, -I- ca 5f^y) , (10) 

oc ((5/^1, - ca 5g^y). (11) 


® This is a generalization of the massive-gravity potential (to 
which it reduces for and a restricted set of ^n) given 

in Ref. [43| . 


^ More general matter couplings not constrained by these require¬ 
ments have been studied in Refs. [^. l44l - [^ . 

^ These correspond to Einstein spaces and, for nonvanishing a, 
solve the field equations only in vacuum. A quartic equation 
determines c = c(/3ti, a). 
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Notice that when a —>■ 0 (or Mpi ^ Mf), the massless 
state aligns with i.e., up to normalization, 

SGfj,i, —>■ + 0{a^). (12) 

Because g^i, is the physical metric, this suggests that 
a —>■ 0 is the general-relativity limit of bigravity@ We 
will see below that the nonlinear field equations indeed 
reduce to Einstein’s equations for a = 0 and that the 
limit is continuous. Thus is close to a GR solution 
for sufficiently small values of a. We therefore identify 
Mpi with the measured physical Planck mass whenever 
a 1, holding it fixed while making Mf smaller. In¬ 
terestingly, in the bimetric setup a large physical Planck 
mass is correlated with the fact that gravity is approx¬ 
imated well by a massless field. In other words, when 
bimetric theory is close to GR, the gravitational force is 
naturally weak. 

The GR limit can be directly realized at the nonlinear 
level [g^, HH . The metric potentials satisfy the identity 

= ^gV5\, (13) 

where V is the potential in the action (El)- For Mf = 0, 
the equation ([6]) gives = 0, an algebraic con¬ 
straint on Then, using the above identity, the 
equation ([5]) becomes 

Gfiuig) + m Vg^n = 2 Tfiv (14) 

Since is conserved, taking the divergence gives 

d^V = 0. (15) 

We see that eq. (0 is the Einstein equation for g^i, 
with cosmological constant m^V. Remarkably, because 
V depends on and all the /3„, this effective cosmologi¬ 
cal constant is generically not simply the vacuum energy 
from matter loops (which is parameterized by /3o)- Even 
in the GR limit, the impact of the spin-2 interactions 
remains and bigravity’s self-acceleration survives. 

It is straightforward to see that, unlike the to —>■ 0 
limit, the a —>■ 0 limit is not affected by the van Dam- 
Veltman-Zakharov (vDVZ) discontinuity [^, [^. The 

cause of this discontinuity is the Bianchi constraint 0 
which constrains the solutions even when m = 0. On 
the contrary, when a —>■ 0, the Bianchi constraint simply 
reduces to eq. 0 and is automatically satisfied. 

The conditions = 0 and = 0 determine 
algebraically in terms of g^i^, generically as = (?g^v 
In the limit Mf = 0, the / sector is infinitely strongly 
coupled0 Due to the nontrivial potential, this causes 


® See Ref. [g^I for an early discussion of such a limit. 

^ Strongly-coupled gravity in the context of GR has been studied, 
for instance, in Refs. IGSHTIH and has been argued to allow for a 
simplified quantum-mechanical treatment. 


the / metric to exactly follow the g metric (both at the 
background and perturbative levels), while the g sector 
remains weakly coupled. 

Strong-coupling scales. — We now argue that at en¬ 
ergy scales relevant to cosmology, this model avoids 
known strong-coupling issues, sometimes contrary to in¬ 
tuition gained from massive gravity. 

There are several strong-coupling scales one might ex¬ 
pect to arise. At an energy scale k, the / sector has an 
effective coupling k/Mf, as can be seen from expanding 
the Einstein-Hilbert action in 6f^^/Mf, just as in GR. 
Then, for small but nonzero a, which is the case of in¬ 
terest here, one might worry that perturbations of 
with momentum k become strongly coupled at low scales 
k ^ Mf. However, we have seen that in the limit of 
infinite strong coupling, M/ = 0, becomes nondy- 
namical and is entirely determined in terms of g^^^ while 
the g^v equation is degenerate with GR and its pertur¬ 
bations remain weakly coupled. Due to the continuity 
of the limit, we expect that, for small enough a, strong¬ 
coupling effects will continue to not affect the g sector, 
even when perturbations of f^i, are strongly coupled at 
relatively small energy scales. In practice, however, since 
the measured value of Mpi is very large, even reasonably 
high values oi Mf can still lead to small a. In cosmo¬ 
logical applications, all observable perturbations satisfy 
k/Mf <C 1 for M/ ^ lOOiJo 10“^^ eV, roughly the 
scale at which linear cosmological perturbation theory 
breaks down at recent times, so that perturbations of 
remain weakly coupled in any case. 

Another potentially-problematic scale is associated 
with the helicity-0 mode of the massive graviton. In mas¬ 
sive gravitY, this mode becomes strongly coupled at the 
scale [zl, 0 

A3= (TO2A/pl)^/^ (16) 

where to is defined to coincide with the Fierz-Pauli mass 
[l| on flat backgrounds. This scale is rather small, 
Ag ~ 10-13 eV ~ (1000 km)-i for m ^ Hq ^ eV, 
and severely restricts the applicability of massive gravity 
[tJ . The same scale also appeys in the decoupling-limit 
analysis of bimetric theory [37| , where to is now the pa¬ 
rameter in front of the potential in the action ([T]). In 
the limit a —> 0, the / sector approaches massive gravity 
and one might worry that the strong-coupling prob¬ 
lem persists or becomes worse with the emergence of an 
even lower scale {m^MfYG^ This is not the case. In 
the bimetric context, the scale defined in eq. 0 is not 
physical, since is degenerate with the /3„. The physi¬ 
cally relevant strong-coupling scale must be defined with 
respect to the bimetric Fierz-Pauli mass [^ . 

TOpp = TO^ ^ c^q;2 T 2c^/32 + Pa) ^ (17) 

which is only defined around proportional backgrounds, 
ffj.v = c^g^iv In the massive-gravity limit, a —>■ oo, the 
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helicity -0 mode is mostly contained in g with a strong¬ 
coupling scale 


A 3 = (18) 

consistent with eq. 0 for appropriately restricted pa¬ 
rameters. However, in the GR limit, a —^ 0, the helicity-0 
mode resides mostly in /, where the strong-coupling scale 
is 


A 3 = ^ ^ > (19) 

which is no longer small. Note that for solutions that 
admit this limit, c becomes independent of a. We can 
also consider the a —?► 0 limit of eq. o, to verify that 
the small part of the helicity -0 mode in g is not strongly 
coupled. 


A3^ ^ . (20) 

This is even higher than A3. Therefore the strong¬ 
coupling issues with the helicity-0 mode are alleviated, 
rather than exacerbated, when a —>■ 0. 

Cosmology. — We now proceed to apply the above argu¬ 
ments to the particular example of a homogeneous and 
isotropic universe. We will take both metrics to be of the 
diagonal FLRW form [l^ - [l^ PI 

g^j.i^dx^dx'' = —dt^ + a'^{t)dijdx^dx^, (21) 

ffj.vdx^dx'' = —X^{t)dt^ + Y'^{t)5ijdx^dx^, (22) 

where we can freely choose the cosmic-time coordinate for 
S/ii' idoo = ~1) because of general covariance. Because 
matter couples minimally to g^i,., this choice is physical, 
and a{t) corresponds to the scale factor inferred from 
observations. We furthermore take the matter source to 
be a perfect fluid, T'); = diag(—/9,p,p,p). The g-metric 
equation ([5|) leads to the Friedmann equation, 

~ 71^2 (^0 + 3/3ij/ -I- 3/322/^ + ! (23) 

iHpj 

where the Hubble rate is defined as H = d/a and the 
ratio of the scale factors is 



The analogous equation for the / metric is 


3A"^ = 




4 + 3% 

y6 y2 


,/53 


+ I 


(24) 


(25) 


® See Ref. and the references therein for other possible metrics 
in bimetric cosmology. 


with K = YjY. The final ingredient is the Bianchi con¬ 
straint 0 , which yields 

(ifX-iF2/)(/3i+2/322/ +432/') =0. (26) 

Taking the first or second term of eq. 0 to vanish se¬ 
lects the so-called dynamical or algebraic branches, re¬ 
spectively. Perturbations in the algebraic branch are 
pathological [ 2 ^, so we will consider the dynamical 
branch in which the /-metric lapse is fixed. 


_ ^ 

H ■ 


(27) 


Inserting this into the equation (I25II transforms it 
into an “alternate” Friedmann equation. 


= w? + 3/32 + 3432 / + 442/'^ • (28) 

We take at least two of the 4n for n > 1 to be nonzero in 
order to ensure the existence of interesting solutions in 
the GR limit a —>■ 0. The solutions to eq. (12^ in the GR 
limit are always on the “finite” branch, i.e., y evolves 
from 0 to a finite late-time value. The perturbations 
on this branch are healthy except for a scalar instability, 
which we discuss below. 

Equation has two features which are useful for our 
purposes. First, in the limit a —>■ 0 it tends to a poly¬ 
nomial constraint that leads to a constant solution for 
2 /, so that the potential term in the Friedmann equation 
(1^51) becomes a cosmological constant. This provides an 
explicit example of the statement above that as a —> 0, 
the theory approaches general relativity with an effective 
cosmological constant (even with 4o = 0). Recall that 
even though the theory approaches GR in this limit, the 
bigravity interactions survive in the form of this constant. 
The other useful feature is that, because eq. does not 
involve p, it can be used to rephrase the potential term 
in eq. (123! ) in terms of the Hubble rate. This will allow us 
to determine the time-dependence of the potential term 
order by order in Q!0 

The effective cosmological constant .— Let us illus¬ 
trate the new viable bimetric cosmologies qualitatively 
by selecting the model with 4o = 43 = /34 = 00 which 
we will refer to as the 4i/32 model. The Friedmann and 
“alternate” Friedmann equations (l23l) and (l28t are 


= iSt + + / 322 /') , 

Mpi 

= + 342 ^ . 


(29) 

(30) 


® One can also combine eqs. an d (l23l to obtain a quartic equa¬ 
tion for y involving p |l34rn. IMI| , but this is more cumbersome 
as it involves higher powers of y than eq. (|2§[) does. 

Since we are interested in finding self-accelerating solutions in 
the absence of vacuum energy, we will set /3o = 0 herein, but 
emphasize that this is not necessary. 
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We can use eq. to eliminate y in eq. (|2^. It is 
instructive to work in the GR limit where eq. (130!) gives 


a^-O 1 Pi 


(31) 


The a —>■ 0 limit is nonsingular only if both /3i and P 2 
are nonzero. Plugging this into eq. (1^ we obtain 


= 


P 


_2^ 2 

Ml , 3/32"" ■ 


The effective cosmological constant is 

A = 


(32) 


(33) 


Late-time acceleration requires /32 < 0. 

When we are not exactly in the GR limit, we should 
consider corrections to eq. 


3iJ^ = 


-b 


Pimp 


Ml , 3 (iL2c,2 _ ^2m2)2 


- 2P2m^) 


2 02 




Ml, 3 P 2 3/?2 
This expansion is valid as long as 
^2 < P^rrP 


Cl(a^ 


(34) 


(35) 



TABLE I. The effective cosmological constant and lowest- 
order corrections (which are time-dependent through p) for 
a variety of two-parameter models. We have chosen solution 
branches which lead to positive Aeft for appropriate signs of 
the jSn, and generally take Pi>G based on viability conditions 
0. The ds, Pi 7^ 0 model does not possess a finite-branch 
solution 11^. 


where $ represents any of the scalar perturbation vari¬ 
ables, N = In a, and we have taken the limit k ^ aH 
where k is the comoving wavenumber. The eigenfrequen- 
cies oj were presented for particular models in Ref. [3^ . 
where it was found that all models with viable back¬ 
grounds have < 0 at early times, revealing a gradient 
instability that only ends at a very low redshift. Using 
the formulation of the linearized equations of motion pre¬ 
sented in Ref. [s^ , we can write the ei genf requencies for 
general /?„ and a in the compact form 41 | 


Rearranging and again keeping terms up to 0{a^), we 
find a standard Friedmann equation with a time-varying 
effective cosmological constant given by 


Aeff 






+ 0(a^). 
(36) 


Because matter is coupled minimally to it will have 
the standard behavior p ~ where w = p/p is 

the equation-of-state parameter, allowing p to stand in 
for time. This captures the first hint of the dynamical 
dark energy that is typical of bigravity [l6l - l^ . 

These results generalize easily to other parameter com¬ 
binations. We list the effective cosmological constant 
up to 0{a^) for all the two-parameter models (setting 
Po = 0) in table Q. We remind the reader that, in order 
for the a —>■ 0 limit to be well-behaved, at least two of 
the Pn parameters (excluding the vacuum energy contri¬ 
bution, Pq) must be nonzero. 


Exorcising the instability. — The stability of cosmo- 
logical perturbations in bigravity was investigated in 
Ref. by determining the full solutions to the lin¬ 
earized Einstein equations in the subhorizon regime. The 
perturbations were shown to obey a WKB solution given 

by 


$ 




( 37 ) 


2 _ (^1 + 4/32^ + 3/33Z/2) y' 

\ k ) ^ iy {Pi + 2P2y + Ps.y'^) 

_ (1 + a'^y'^) {Pi - Psy'^) 

3a2y3p(i + 

where p = p/mpMl, and primes denote d/d In a. 

We apply this to the P 1 P 2 model. Assuming a universe 
dominated by dust (w = 0), crosses zero whei0 

180^/32 {a^Pl + ^Pl) y^ + ^a^Pi {a^Pl + 10 / 3 |) y^ 

+A8a^PlP2y'^ + 6P2 {‘^a^Pf - Pi) y^ - &PlP2y - Pl = t). 

(39) 

Solving this for y, we can then use eq. to determine 
the value of Hubble rate at the transition era., before 
which the gradient instability is present and after which 
it vanishes. While this solution is too complicated to 
write down explicitly, in the limit a —7 0 the leading- 


We have used e ns. ||29| '| and 1301 1 and their derivatives to solve for 
y' and p in eq. J38f) in terms of and y l3ll . Note that = 0 
does not imply strong coupling because, while the gradient terms 
vanish, the kinetic terms remain nonzero. 



























6 


Era of transition to stability 


a 

Mf 

BBN 

IQ-i® eV 

10-1" 

100 GeV 

A3 = (m^Afpi/a) ^ 

10“® eV 

10-^1 

10“® eV 

GUT-scale inflation 

GeV 

10-55 

10-2" eV 

Mpi 

10^® GeV 

10-61 

10-^® eV 


TABLE II. The values of a and Mf for a few choices of the 
era at which perturbations become stable. 


order term is remarkably simpleF^ 

Hl = ±^ + 0{a% (40) 

where iT* is defined as the Hubble rate at the time when 
= 0, i.e., after which the gradient instability is absent. 
We pick the negative branch of eq. for physical rea¬ 
sons, i.e., so that > 0 given that /32 < 0. We have 
checked explicitly that by solving for y with this value of 
H and plugging it into all terms up to 0{a^) vanish. 

Interestingly, eq. o is the same as the condition (1351) 
for the small-a expansion of the background solution to 
be valid. Therefore, simply by pushing the instability 
back to early times, one gets late-time bimetric dynam¬ 
ics that can be described as perturbative corrections to 
GR, except for the effective cosmological constant which 
remains nonperturbative. This is nontrivial; while we ex¬ 
pect everything to reduce to GR at late times when we 
can expand in aH /there could in principle have 
been earlier times during which perturbations were stable 
but still fundamentally different than in GR. 

We can rewrite eq. in more physical terms as 


Hi 


3^3 



(41) 


where is the far-future value of H and should be 
comparable to the present Hubble rate, Hq. For |/?i| 

1/321, this implies simply 




a 


(42) 


We see that as we approach the GR limit, the smaller 
one takes the /-metric Planck mass, the earlier in time 
bigravity’s gradient instability is cured. Our goal is to 
make this era so early as to be effectively unobservable. 
One has a variety of choices for the scale where the insta¬ 
bility sets in; the values of a and My for various choices 
are summarized in table |l3. 


While eq. Hdi only holds exactly in the presence of dust, w = 0, 
for other reasonable equations of state, such as radiation {w = 
1/3), it will only be modified by an 0(1) factor. Since we will be 
using this analysis only to make order-of-magnitude estimates, 
the exact factors are unimportant. 


A natural requirement would be to push the instability 
outside the range of the effective field theory, i.e., above 
either the cut-off scale where new physics must enter, or 
the strong-coupling scale where tree-level unitary breaks 
down[3 The cut-off scale in massive and bimetric gravity 
is not known. The strong-coupling scale, to the extent it 
is understood, was discussed above. Here we focus on ob¬ 
servational constraints. It is natural to demand that the 
instability lie beyond some important cosmic era which 
we can indirectly probe, such as big-bang nucleosynthesis 
(BBN) or inflation. Both of these possibilities are then 
likely to be observationally safe as long as the Universe is 
decelerating (e.g., is radiation-dominated) after inflation, 
because the instability is only a problem for subhorizon 
modes with large k/aH, and during a decelerating epoch 
modes with fixed comoving wavelength always become 
smaller with respect to the horizon. Consider, as an ex¬ 
ample, that the transition to stability occurs between 
inflation and BBN. During that period, modes will grow 
rapidly on small scales, but those will be far, far smaller 
than the modes relevant for the cosmic microwave back¬ 
ground or large-scale structure. One might worry that in¬ 
flation’s ability to set initial conditions is spoiled in this 
scenario (assuming that the linear theory is even valid 
during inflation, which is not guaranteed due to the argu¬ 
ments above). However, the instability should be absent 
during inflation; notice from eq. (fH) that generically 
becomes large and positive for w close to — iE3 Therefore 
the instability would not affect the generation of primor¬ 
dial perturbations during inflation. If the instability later 
appears with the onset of radiation domination, it would 
only affect small scales which are irrelevant for present- 
day cosmology. 

If the instability ends at the time of BBN, Mf can be 
as high as about 100 GeV, far larger than the wavenum¬ 
bers probed by cosmological observations. We remind 
the reader that for such a “large” Mf, perturbations in 
the Einstein-Hilbert term for remain weakly-coupled 
for all observationally-relevant k. 

While analytic results like eq. 0 cannot be obtained 
for most of the other two-parameter models, we have 
checked that in each case the relevant behavior, ^ 
H\ja^ holds 13 The values given in table [H are therefore 
fairly model-independent. 

The other pathology that is typical of massive and bi¬ 
metric gravity, the Higuchi ghost, is not present in these 


These two are not always the same, and may not be in massive 
and bimetric gravity [70 .17^ . 

This depends on the exact parameters and the evolution of y. 
Background viability requires /3i > 0 [l9ll . so as long as ^3 < 0, 
at the very least the last term in eq. |^) is large and manifestly 
positive. 

Specifically, this holds in the models with /3i 7 ^ 0. The gradient 
instability is absent from the I32l3s and / 32/34 models at early 
times [32| . These were shown in Ref. [l3| to have problematic 
background behavior at early times, but these again can be made 
unobservably early in the GR limit. 
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models. There is a simple condition for the absence of 
this ghost, dp/Ay < 0 (H, [s^ (see also Refs. [Hill). 
Because for normal matter p is always decreasing with 
time, this amounts to demanding that y be increasing. 
In the “finite-branch” solutions which we are considering, 
y evolves monotonically from 0 at early times to a fixed 
positive value at late times, and so the Higuchi bound is 
always satisfied [4lj| . 

Parameter rescalings. — We have presented a physi¬ 
cally well-motivated region of bimetric parameter space, 
near the GR limit, in which observable cosmological per¬ 
turbations are stable and yet self-acceleration remains. 
One is naturally led to ask how this has been missed by 
the many previous studies of bimetric cosmology. The is¬ 
sue lies in a rescaling which leaves the action o invariant 

[Mill, 

fnu —t /?„ —>■ -^Pn, Mf —>■ flMf, (43) 

and hence gives rise to a redundant parameter. It has 
become common to let a play this role and perform the 
rescaling O = l/a such that a is set to unity. While 
our results do not invalidate this rescaling, they do show 
that it picks out a particular region of parameter space 
which may not capture all physically-meaningful situa¬ 
tions. In particular, the a —>■ 0 limit, in which the theory 
approaches GR—the behavior at the heart of our remov¬ 
ing the gradient instability—would look extremely odd 
after this rescaling: the /3„ would not only be very large, 
but each /3„+i would be parametrically larger than /3„ll£| 
Therefore, studies which set a to unity could in principle 
have found the GR-like solutions which we study here, 
but only by looking at what would have appeared to be 
a highly unnatural and tuned set of parameters, even 
though they have a simple and sensible physical explana¬ 
tion. Without performing this rescaling, we can simply 
take the nonzero /3„ to be 0{1) and consider that we are 
in the small-Mj regime. 

It is clear that in phenomenological studies of bigrav¬ 
ity, a must not automatically be set to unity. When 
working with a two-/3„ model, perhaps a more sensible 
rescaling would be one such that the two /3„ are equal 
to each other (up to a possible sign). They can further 
be absorbed into In this case, the free parameters 
are effectively the spin-2 interaction scale, m^, and the 
/-metric Planck mass, Mf. Their effects decouple nicely: 
Mf controls the earliness of the instability, while m sets 
the acceleration scale. Alternatively, one may consider 
that the rescaling (l4^ simply tells us that rather differ¬ 
ent regions of parameter space happen to have the same 
solutions, and therefore not perform any rescaling a pri¬ 
ori at all. 

Summary and discussion. — We have shown that a 


We can recast this as a large m^, but there would remain a 
specific tuning among the jSn of the form finlPn+l ~ e, where e 
is the value of a before the rescaling. 


well-motivated but heretofore underexplored region of pa¬ 
rameter space in bimetric gravity can lead to cosmolog¬ 
ical solutions which are observationally viable and close 
to general relativity, with an effective cosmological con¬ 
stant that is set by the spin-2 interaction scale m. In 
this limit, obtained by taking a small /-metric Planck 
mass, the gradient instability that seems to generically 
plague bimetric models at late times is relegated to the 
very early Universe, where it can be either made unob¬ 
servable or pushed outside the regime of validity of the 
effective theory. This instability had been considered in 
previous work to make bimetric cosmologies nonpredic- 
tive even at late times. Furthermore, in this limit the 
theory avoids the usual low-scale strong-coupling issue 
that affects the helicity-0 sector in the massive-gravity 
limit. 

What is encouraging is that the one property of bi¬ 
gravity which survives in the small-a limit is its cosmo- 
logically most useful feature, the technically-natural dark 
energy scale. In other words, the effective cosmological 
constant of bigravity in a region close to GR is not just 
the vacuum-energy contribution and can give rise to self¬ 
acceleration in its absence. 

The model we have presented is expected to be ex¬ 
tremely close to GR at all but very high energy scales. 
In particular the Newtonian limit is well-behaved; unlike 
—>■ 0, which suffers from the vDVZ discontinuity, the 
GR limit a —>■ 0 is completely smooth because all the 
helicity states of the massive spin-2 mode decouple from 
matter. Note also that massive gravity does not possess 
such a continuous GR limit. 

It is worth emphasizing that the a —^ 0 limit brings 
bimetric theory arbitrarily close to GR even for a large 
value of the spin-2 mass scale, m ^ Hq. The presence 
of heavy spin-2 fields in the Universe is therefore not 
excluded as long as their self-interaction scale (set by 
Mf) is sufficiently small compared to Mp\. In this case, 
however, the /3„ parameters need to be highly tuned for 
the effective cosmological constant small enough to be 
compatible with observations 0 Note however tha.t. since 
the Pn are protected against loop corrections 0, , 

this tuning does not violate technical naturalness. 

Finally we comment on the potential observable signa¬ 
tures of this theory. While at low energies, corresponding 
to recent cosmological epochs, this limit of bigravity is ex¬ 
tremely close to GR, there may be observable effects at 
early times when the effects of strong coupling become 
important. In this case, given hy H > i7*, the small-a 
approximation breaks down and modified-gravity effects 
must be taken into account. This may be particularly 
important for inflation, which will see such effects unless 
Mf is extraordinarily small. A better understanding of 


Indeed, without this tuning of the fin, the interaction term would 
lead to acceleration at an unacceptabl y e arly epoch. This sce¬ 
nario is related to the findings of Ref. 1361 , where it was shown 
that the instability becomes negligible for large values of m. 




strong coupling in the /^j/ sector will therefore point the 
way towards tests of this important region of bimetric 
parameter space, since at this point it is not clear how 
to perform computations in the strong-coupling regime. 
There may also be effects related to the Vainshtein mech¬ 
anism [80, [MI- We conclude that the closeness of this 
theory to GR is both a blessing and a curse: while it is 
behind the exorcism of the gradient instability and brings 
the theory in excellent agreement with experiments, it 
presents a serious observational challenge if it is to be 
compared against GR. It is nevertheless encouraging that 
this “GR-adjacent” bigravity naturally explains cosmic 
acceleration while avoiding the instabilities that plague 
other bimetric models, and therefore merits serious con¬ 
sideration. 
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